A momentum-space representation of Green's functions 
with modified dispersion on ultra-static space-time. 
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We consider the Green's functions associated to a scalar field propagating on a curved, ultra- 
static background, in the presence of modified dispersion relations. The usual proper-time deWitt- 
Schwinger procedure to obtain a series representation of the Green's functions is doomed to failure, 
because of higher order spatial derivatives in the Klein-Gordon operator. We show how to overcome 
this difficulty by considering a preferred frame, associated to a unit time-like vector. With respect 
to this frame, we can express the Green's functions as an integral over all frequencies of a space- 
dependent function. The latter can be expanded in momentum space, as a series with geometric 
coefficients similar to the deWitt-Schwinger's ones. By integrating over all frequencies, we finally 
find the expansion of the Green's function up to four derivatives of the metric tensor. The relation 
with the proper-time formalism is also discussed. 

PACS numbers: 04.60.m 
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I. INTRODUCTION 

Modified dispersion relations (MDRs) have recently at- 
tracted a flurry activity in various high-energy physics 
models. The common motivation is that MDRs can 
be used as a phenomenological approach to investigate 
physics at the Planck scale, where General Relativity is 
no longer reliable. In general, the lack of a complete the- 
ory of quantum gravity, leads to consider trans-Planckian 
effects as perturbative above a certain energy scale. The 
first application is certainly in cosmology. Several au- 
thors believe that trans-Planckian effects affected the 
early-stage evolution of the Universe, and that they left 
some observable fingerprints, e.g. in the CMB inhomo- 
geneities (see, for example, [l|). A similar situation oc- 
curs in black hole physics, where the quantum thermal 
emission discovered by Hawking is related to modes 
of arbitrarily large frequency near the horizon. In this 
case, it has been proved that the spectrum emitted at 
infinite distance from the hole is only marginally affected 
by MDRs 3]. 

Some fundamental theories, such as String Theory, 
Loop Quantum Gravity, and Double Special Relativity, 
predict MDRs (see the review [4*1 for references). How- 
ever, little is known about the modifications that MDRs 
generate in the formalism of quantum field theory on 
curved space. For example, the renormalization of the 
stress tensor is crucial to evaluate the back-reaction in the 
semi-classical theory [5] . Such a quantity, in the presence 
of MDRs, has been recently obtained in the context of 
cosmology, by the authors of [^] . In this work, it is shown 
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that the renormalization procedure leads to a rescaling 
of the bare Newton's constant and cosmological constant. 
However, in the case of non-homogeneous backgrounds, 
the renormalization of the stress tensor, in the presence of 
MDRs, appears much more difficult, as the Klein-Gordon 
operator, unlike the cosmological case, now contains spa- 
tial derivatives of (at least) fourth order (see for example 
Q). In fact, the usual deWitt-Schwinger representation 
of the Green's functions, which is the starting point of 
the point-splitting technique, [l]-[l0| does not work in 
this case. 

In this paper, we begin to consider the problem of 
MDRs on non-homogeneous manifolds by first looking at 
the ultra-static case. Our method, which should be appli- 
cable also to stationary metrics, relies upon the existence 
of a preferred frame. It is known that, even though MDRs 
break the local Lorentz invariance, general covariance can 
be preserved by introducing a preferred frame through a 
dynamical unit time- like vector field [iH-[l2|- With 
the help of the latter, we can foliate the manifold into 
space-like surfaces. Furthermore, if the space-time is also 
ultra-static, the Green's functions can be written as an 
integral over all frequencies of a function, which is inde- 
pendent of the time associated to the observer co-moving 
with . This function satisfies an equation, which can be 
solved in momentum space by applying the well-known 
Bunch and Parker procedure [13l |. Thus, we can find a 
momentum-space representation of the time-independent 
part of the Green's functions. 

The plan of this work is the following: in the next sec- 
tion, we introduce the modified dispersion relation which 
will be used in this work, and, in section 3, we calculate 
the two-point function in flat space. In section 4 we intro- 
duce the deWitt-Schwinger analysis, suitably adapted to 
our case. In section 5 we present the Bunch and Parker 
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method, and we find the expansion of the Green's func- 
tion in momentum space up to four derivatives of the 
metric. Finally, we conclude with some remarks and fur- 
ther conjectures. Throughout this paper, we use the sig- 
nature (— , -|-, -f • • • ), and set = c = 1 



II. MODIFIED DISPERSION RELATIONS 

In Minkowski space-time the dispersion relation for a 
scalar field of mass m can be found by inserting the 
ansatz (f> ~ exp{—ikot + ik- x) in the Klein-Gordon equa- 
tion 

(-V^+to2)0^O, (1) 

which, together with the identity ujk — |fco|, yields 

u;l^\k\^+m^. (2) 

A general dispersion relation can be written as 

u;l^\k\'+H\k\)+m' , (3) 

where is a scalar function of the modulus of the 

wave-vector k. If depends on the square of the modu- 
lus, rotation invariance is preserved. If it is also analytic, 
then it can be expanded, and, to leading order, the MDR 
reads 



Ifcl 



(4) 



where e is a cut-off parameter that sets the lowest value 
of k at which corrections to the dispersion relation are ig- 
nored. Also, the sign of indicates wether the modes are 
sub- luminal (e^ < 0) or super- luminal (e^ > 0). These 
kinds of MDRs were considered in Cosmology, but also in 
the context of the analogue models of gravity constructed 
with superfluids [1]. It is clear that MDRs, such as 
break Lorentz invariance. However, general covariance 
can be preserved if the preferred frame is associated to a 
dynamical quantity. This is the route followed by Jacob- 
son et al. [ll|: the preferred frame is determined by a 
unit timelike vector field which enters quadratically 
the action. The latter has the form 



S ■ 



aiR - biF^'^Fab + HgabU^u'' + 1)) 



(5) 

where oi and bi are constant parameter. Fab — 2V[QMf,], 
and A is a Lagrange multiplier that ensures to be a 
unit time-like vector. It is then possible to construct a 
massless scalar field Lagrangian, which preserves general 
covariance, and reads 



r = 



(6) 



where e sets the scale at which Lorentz invariance breaks. 
The operator is the covariant spatial Laplacian de- 
fined as 



V^^ = -q--Va{qJ'VbV) , 



(7) 



where qab is the induced metric on the spatial sections, 
defined by 



qab ^ gab + UaUb 



(8) 



The inclusion of the vector field m° as a dynamical vari- 
able leads to a theory similar to the Maxwell-Einstein ac- 
tion. Therefore, one expects an Abelian conserved charge 
associated to However, it is possible to show that if 
one assumes a FLRW flat metric, the equations of motion 
are trivially solved, i.e. Fab = 0. On the contrary, sym- 
metries do not allow to choose Fab = 0, together with 
a Schwarzschild metric. In fact, the conserved charge 
associated to Fab does not vanish, leading to a Reissner- 
Nordstrom form for the metric [HJ. Therefore, when 
MDRs associated to dynamical vector fields are consid- 
ered, the modes propagate on a modified metric. In this 
case, one can argue that the modifications introduced 
by can be somehow set to be very small, so that the 
modes with a MDR can still be seen as propagating on 
an unchanged background. In any case, in order to eval- 
uate the effects of MDRs at the level of back-reaction, 
we first need to compute the renormalized Green's func- 
tions. In order to do so, we begin by evaluating the 
deWitt-Schwinger expansion in section 4 and 5, while in 
the next section we look at the simple flat space case. 



III. MDRS IN FLAT SPACE 

Before embarking upon the study of the curved back- 
ground case, let us have a look flrst to what happens in 
flat space. For simplicity, we consider the 2-dimensional 
Minkowski metric 



ds^ ^ -df + dx^ 



(9) 



where the Klein-Gordon equation, obtained from the La- 
grangian reads 



(10) 



Despite the modification, the standard function 
exp(— iwt -|- ipx) is a basis for the (physical) solution 
space, provided uS^ = + e^p*. This might be confusing, 
since we expect four independent solutions. Indeed, this 
is the case but two out of the four solutions have imag- 
inary momenta. To show this, it is sufficient to replace 
the ansatz (f)(x,t) ~ (^(a;) exp(— iwpt) into Eq. (jlOp . and 
find the fourth-order differential equation 



e^d^(p{x) - dlip{x) - uj^Lp{x) = 



(11) 



The latter has four independent solutions of the form 

ViA^) = Ci,2 , (^3,4(a;) = 6-3,4 , (12) 

where Ci,2.3,4 are normalization constants, and 



1 ± 



V2e 



(13) 
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We note that fl^ Ms a pure imaginary number and that, 
for e — > 0, it tends to the usual factor iLUp. On the con- 
trary, f2(+^ is real and it diverges when e — )■ 0. By solving 
the above identity for LUp, one finds that 



(±)^4 



(14) 



which is nothing but the modified dispersion relation, 
provided we identify fi*^"' with ip. It then follows that 
= p and the corresponding functions (/?3^4 span an 
unphysical sector of the solution space [T7| . 

In momentum space, the massless propagator corre- 
sponding to the modified dispersion reads 



G{p) 



1 



1 



UJt, 



Pi 



(15) 



To compute the two-points function, it is convenient 
to find first the Wightman functions and G~, by 
choosing the appropriate contour C in the 2-dimensional 
Fourier transform 181 



G(a;^,a;'^) 



dpcidp e 



Po 



(16) 



where Ax^ 



The pole structure of the above 



integral is formally the same as in the relativistic case: by 
adding a small imaginary part to the denominator, the 
poles are located at pq = ±ujp =p iS, with S << 1. The 



functions iG~^ and iG~ can be found directly by choosing 
C as circles around each pole (clockwise around —tUp and 
counter-clockwise around +ujp, see [5]). The result is 



G*(a;^,2;'^) 



2tt 



cos(pAx) ^(zpiAt^/^^H^v)^^ 

(17) 

where A is the usual IR cut-off, and Ax = x' — x, 
At ~ t' — t. The breaking of Lorentz invariance is now 
manifest, since Ax and At do not have the same p pre- 
factor. This integral cannot be solved exactly, hence 
some approximations need to be done. First, we note 
that the ratio between the integrand in Eq. p7|) and the 
integrand of the relativistic propagator 



1 



(x^x"') = - 

ZTT 



+ 00 



cos(pAx) 



P 



'S^'P^'^dp, (18) 



tends to one as p ^ 0+ . It is then reasonable to assume 
that the IR behavior of the modified Green's function is 
the same as the relativistic one. Therefore, we do not 
worry too much about the IR divergence and we look at 
the large p regime, by approximating the integral (fTT]) as 



d 1 /■+°° 

—— (G+ + G~) = / dpcos(pAx) sm(ep'^At) 

dAt 2tt Jq 



which yields 



(19) 



=G+ + G- 




f Aa 



V4eAi 



/ Aa 



\UAt 
I 



Ax 
"27 



Vv27reAi 



Ax 



V27reAt, 



, (20) 



where C and S are Frenel's integrals [19| . It is interest- 
ing to note that in the coincidence limit = x'^, this 
expression is actually finite. However, we expect diver- 
gences in the coincidence limit of the energy-momentum 
tensor, as the latter is calculated trough derivatives of 
the two-point function with respect to At and Ax. 



IV. THE DEWITT-SCHWINGER EXPANSION 

In the previous Section we found an expression for the 
two-point function at large momenta and in flat space. 
In the relativistic case, the generalization of such ex- 
pressions on curved space proceeds through the deWitt- 
Schwinger construction [8]- [10], which works directly in 
coordinate space. Here we briefly outline the procedure, 
adapted to the modified dispersion relation (|4]), with 

> 0. We consider a. n + 1-dimensional globally hyper- 
bolic manifold, such that it can be foliated into space-like 
surfaces of constant r [l3 | . The parameter r can be used 



to define the unit time-like tangent vector Up = — 9^t, 
with respect to some coordinate system x^. Therefore, 
T assumes the role of the time relative to a free-falling 
observer moving with velocity u'^. Thus, we can write 
the metric as 



ds^ — gpiydx^dx^ = — dr^ + q^^dx'^dx'^ 



(21) 



where dr — Updx^. In this paper, we consider the case 
when q'"' = g^" + u'^u'^ does not depend on r, i.e. when 
the metric is ultra-static The form of this metric guaran- 
tees that 



detg^i, = detqij 



(22) 



where Latin indices, here and in the following, label spa- 
tial coordinates only. Thus, the d'Alambertian operator 
acting on a scalar can be decomposed as 

VV = V20-0, (23) 

where the dot indicates derivative with respect to r. 
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Modified dispersion relations (|4]) naturally appears if 
one considers the Klein- Gordon equation 



where we changed s is, and set a 



(24) 



for which the associate Green's functions satisfies the 
equation 

{V^^m^ -e'V^)g{x,x') ^ -g-^'H^^+^\x^x') . (25) 

Here, and from now on, we denote g{q) = det |(7(g)|. Be- 
cause of the MDR, the function Q is no longer Lorentz 
invariant, but, since the modification is a quartic spa- 
tial operator, the 0{n) rotational invariance holds. It is 
therefore convenient to express the Green's functions as 

G{x,x') = j '^e'^^^-^'^G{x\x'\uj) , (26) 

where G{x^ , x'^ , uj) depends on spatial coordinates and lj 
only. With the help of the identities ([22]) and ([23]), one 
can easily show that G must satisfy the equation 

(V^ - + - £2 v4) G{x' , x'^ , ^) ^ 

= _q-l/2j(n)(2.J ^ (27) 

provided the metric ((2T|) is assumed. 

This equation describes the Green's functions associ- 
ated to a scalar field, which propagates on the spatial 
sections of the metric (|211) . as seen by a free- falling ob- 
server co-moving with u'^. From the mathematical point 
of view, we can construct the heat kernel functional 
H{s,Ax^), formally defined as 



pec 

G{Ax^,Lj) = / H{s,Ax^)ds . 
Jo 



Then, H must satisfy the heat equation 
OS 



(28) 



(29) 



where Z is an elliptic operator. In our case, Z — — 
TO^ + - e'^W*. 

In the relativistic case, e = 0, one can show that the 
heat kernel has the form [T3| 

H{s, Ax\e = 0) - (47rs)-"/2 ^-i^x^ /is+s(u.-'-n{') 

X ^ ai{x^,x'^)s^ , (30) 

(=0 

where n is the number of space-like dimensions and the 
a' s are geometric coefficients built on the Riemann tensor 
associated to q^v, which can be found by a recursion 
procedure (see @). Finally, by integrating over w, one 
find the full Green's functions expansion 

/•oc 

g{x, x') = i / ids(4^zs)-("+i)/2 ™' 



Y,ai{x^,x'^){zsy , (31) 

1=0 



A crucial remark is in order here. Despite the above 
formula looks identical to the usual deWitt-Schwinger 
expansion, there are major differences. First, the a's 
depend on geometrical quantities built on g^i,, not 5^^. 
This is due to the fact that the metric ((2T|) is not the 
most general one, as the space-time under consideration 
is static with respect to the free-falling time r. As a con- 
sequence, the frequency w is not a generic one, but it is 
the one measured by the free-falling observer only. 

To obtain the deWitt-Schwinger expansion in terms of 
geometrical quantities related to g^^^, some extra work is 
required. As an example, consider the expansion up the 
the second order only. The only non-trivial coefficient 
reads [lol | 
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(32) 



where R is the n-dimensional Ricci scalar built on q^^ . 
By contracting the Gauss-Codacci identity [3| 

R'^p^s = ^ W?" 7? 5 - K%K^s + K^sKp-t , (33) 

where K^p — q^^q^pS^s Uj is the extrinsic curvature of 
the spatial sections orthogonal to u^, and R'^^^p is the 

Riemann tensor constructed with 5^1^ , we can relate R to 
the curvature of the full space-time, as 



R = R + 2R^fiU°'u'^ ~ {K°'J' + K°"'K^ 



Q/3 j 



(34) 



Thus, Eq. (I31|) contains extra terms with respect to the 
usual expression [9] 

/"OO 

g{x,x') = I / ids(47rzs)~("+i)/2e-'^/2»s-».m^ 



J2aiix,x')iisy . (35) 



1=0 

We now turn to the dispersive case. When e ^ 0, the 
ansatz (1501) simply does not work, and the reason will 
be clear at the end of next section (see also [l^l)- As 
our 71-dimensional problem can be seen as the Euclidean 
continuation of a n-dimensional theory with a term, 
one can find the heat kernel by the method presented in 
p^ . However, as we will see in the next section, there is 
an alternative route, which simplifies the calculations and 
keep track closely of the physics behind the mathematical 
structure. 



V. MOMENTUM SPACE REPRESENTATION 

In this section, we find a momentum-space represen- 
tation of the Green's functions, which, in principle, can 
work for any analytic dispersion relation J-'{k^). For the 
moment, we restrict ourself to the quartic case and 
discuss further generalizations in the last section. The 
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method that we are going to use is due to Bunch and 
Parker [l^, and it essentially makes use of a local 
Taylor's expansion of the metric tensor expressed in Rie- 
mann normal coordinates 0, [l^ . In turn, this leads to 
a similar expansion of the various differential operators 
in terms of ordinary partial derivatives and geometrical 
coefficients. To begin with, we define the function G such 
that 

G{x, x',Lo)^ q-^'^{x)G{x, x\ uj)q~^'^{x') , (36) 
and write 

q-i/2(a;)5(x - x') = q-^'^(x)5{x - x')q-^'^{x') . (37) 

These function behaves as a bi-scalars at x and x' . For 
notational convenience, we drop the indices, and x and 
x' are considered as separate point on the same spatial 
slice T — const. Next, we introduce the Riemann nor- 
mal coordinates y with origin at x' . Thus, q{x') — 1 
and the Green's function G satisfies (from now on, the 
dependence of G from lo will be understood) 

+m^-V^+ e2v^)(g-i/4G(y)) ^ S{y) . (38) 



(43) 



where 



^ 30 60 ''^^ '^ 60 i^pqij 

+ 40^:^^ + 40 ^^^'P^ ■ 
The expansion above is also obtained by using the fact 

that G depends on y^, being rotationally invariant. To 
find the expansion for the quartic operator, we can pro- 
ceed by iteration, but first it is convenient to move into 
momentum space. Hence, we define the local Fourier 
transform of G as 



G{y) = 



(27r)" 



G{k) 



(44) 



Now, let 



^ = gi/4v2((7-i/4G(y)) . 



(45) 



In momentum space, this corresponds to 



In a neighborhood of x' , we can expand the induced met- 
ric as uM 



Qmn — ^mn 
1 



1 



:Rn 



-.hy^y^ - lRmanb;py°-y''y^ 

6 



20 



Rjnanb-.pq 



y'^y'y^y'^ , (39) 



from which it follows that 

9 - 1 - l.Raby-'y' - \Rab,cy''y'y' 
3 6 



— RabRcd - -^Rpab'^R^cdq ~ l^Rab-cd^ y°'l/'v''v'^ ■ 

(40) 



All coefficients are evaluated at x' (i.e. at y = 0) and 
contain up to four derivatives of the metric qij. If we 
write 



V; = -k^G+-RG+-R.ikWG-2HDG + 
6 3 

- AHijk'kW^G , (46) 



where we defined the operator D such that 



In coordinate space, V^{q ^^'^G{y)) = q ^/''^ and 

5l/4v4(^-l/4(5) ^ 5l/4v2(^-l/4^) _ (48) 

As also q~^^^^p is rotationally invariant, all we need to 
do is to insert the expansion ((46|) into itself, to find that 
the Fourier transform of q^^^V^{q^^/^G) is given by 



gi/4v2(q-i/4G(2;)) = q^3d,d,G + d,q'^d,G 
"^g'^a,(lng)a,(lng) + 

-^q'^d,d,{\nq) + -^d,q'^d,{\nq)\G , 



(41) 



1,2 1 1 

^ k'^G RG - -Rjk^G + 2k'^DG) + —R^G 



36 

2H{G + 2PDG) + ^H,jk^ki{k^D'^G + DG) 



(49) 



we can evaluate all the coefficients and find, up to fourth 
order. 



?V4v2(q-l/4G(y)) 



5''d,djG + ^RG + ^kjfG + 



H^,yYG , 



(42) 



where H = 5^^ Hij . Note that this iterative procedure 
can be used to find the expansion, in momentum space, 
of for any integer p. 

With these elements, we can expand Eq. in mo- 
mentum space up to four derivatives of the metric as 



6 



6 



-R,jk^{DG + e^G + 2e^k^DG) ' ^ 



36 



2i? (DG + e^G + 2e^k^DG) + AH^k^y {D^G + 2€^k^D^G + 2e^i:>G) = 1 



At the zeroth order, this equation yields 
Go 



fc2 + e^k^ + — Lo'^ 



(50) 



while the following orders can be found by recursion, 
yielding 



1 



(51) 
(52) 



G2 = -i?(l + 2e^A:^)G^ , 

G3 = iAjFGo(G + 2A;2i^Go) , 

G4 - ^R%l + 2e^efGl-^^R^Gl + 
+ 2i7Go(I?Go + 2e2A:2£)G<g + g2(5^)_^ 

+ 4H,jk'k^Go{D^Go + 2e^k^D^Go + 2e^DGo) . 

(53) 

Therefore, as G = Go + G2 + . . ., we finally have 

G = Go - ii?i?Go - -^A.a-'.DGo + 
6 12 

^R" - ^H^ D^Go - ^H^.d'dWGo , (54) 

where, tilded derivatives are with respect to fc% and 
where many cancellations occur by using the identities 

i^Go = -{l + 2e^P)Gl , (55) 
(1 + 2e2fc2)D2G2 = -4e2DG^ - D^Gq , (56) 



are built on qtj . These coefficients are formally identical 
to the ones found by Bunch and Parker in [l3|, but, in 
that case, they were built on the full metric (7^1/. It is 
interesting to note that, when, = we have 



d 
9m 2 



(61) 



and the above expression becomes identical to the Bunch 
and Parker one, except, again, for the geometrical coeffi- 
cients. As mentioned above, this is due to the fact that 
we are working with the metric (|2ip . which is not the 
most general one. 

At uj fixed (i.e. on a given spatial slice), we define 



2_]_,2l,4 I „,2 ,2 



io 



(62) 



and replace into ([54]) . By swapping integrals, we find 



l + {is){l~2e'd')h + 



G{y) = i / dse-*''(™ ) 

+ 2{isfe^h + {isf{l~2e'd^ff2]h{y.s), (63) 
where the partial derivatives are with respect to y, and 
I^iv^^)- |^e-('=^+^^'=^)+^'=-V (64) 

This integral can be evaluated in form a sum of Hermite 
polynomials 7i[;] (see the Appendix), namely 



Goi^^Go = ^D^Gl - '^e^Gl . (57) 



The expansion ((54l) is our main result. We conclude 
this section, by connecting the above expansion with the 
proper time formalism of deWitt-Schwinger. By using 
Eqs. by setting w = fc", and by integrating 

by parts, we find the full (n + l)-dimensional Green's 
function, expressed as 



where 



(58) 



/i = ^R;jy' - lH,,fy' , (59) 



(60) 



_ V- f — 

\n/2 \| I 1R„ 



(4is7r)"/2 ^ A! 

^ ' A=0 



16s 



n 



[4A] 



/4is 



(65) 

It is clear that the deWitt-Schwinger expansion in proper 
time (j63| becomes a very complicate sum of Hermite 
polynomials and their derivatives, which does not appear 
to converge to any known function. In the Lorentz invari- 
ant case instead, the sum is trivial and the integral over s 
converges to Hankel functions of second kind [1, [ol] . It is 
now clear why the heat kernel ansatz l|30p cannot work. 



VI. DISCUSSION 

In this work we present our first results concerning quan- 
tum field theory on curved backgrounds, with modified 
dispersion relations. In particular, we carefully analyze 
the Klein-Gordon equation associated to a scalar field 
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propagating on a ultra-static space-time, as a first step 
towards the physically relevant case of stationary met- 
rics. 

In the simplest case, superluminal propagation of high 
frequency modes can be achieved by adding a quartic 
spatial derivative to the Klein-Gordon operator. An im- 
portant side-effect of these higher order operators is the 
breaking of Lorentz invariance. On physical grounds, 
this leads to important modifications to the Green's func- 
tions, even in flat space. As an example, we considered 
the 2-dimensional flat space, and we proved that the two- 
point function becomes finite in the coincidence limit. 

Motivated by this, and also by the results achieved 
in the context of Cosmology, we decided to consider the 
problem in a non-homogeneous curved space. The prob- 
lem becomes much more difficult, as in the Klein-Gordon 
equation we now have higher-derivative spatial operators, 
in contrast to the cosmological case, which prevent the 
usual deWitt-Schwinger expansion of the Green's func- 
tions from working properly. We turned around this 
problem by assuming the existence of a preferred frame, 
encoded by the unit time- like vector u^, associated to the 
free-falling observer. In this way, despite the breaking 
of local Lorentz invariance, we still have general covari- 
ance and, above all, rotational invariance over slices of 
constant time. In the case of a ultra-static space-time, 
this fact allows for a dimensional reduction of the Klein- 
Gordon equation, so that one can consider unambigu- 
ously the Green's functions at a fixed frequency, as mea- 
sured by the free-falling observer. As a result, one deals 
with elliptic operators at constant time. Despite this sim- 
plification, the usual heat-kernel ansatz does not work yet 
because of the higher derivative operators, responsible for 
the modified dispersion relation. Therefore, it turned out 
to be much easier to work in momentum-space, follow- 
ing the method of Bunch and Parker. Thus, we obtained 
an expansion of the Green's functions up to four deriva- 
tives of the metric in momentum space. We also showed 
that the proper-time formulation of our expansion leads 
to a very complicate expression, which might render very 
difficult to renormalize the theory in coordinate space. 

We wish to conclude with a remark. In this work, 
we displayed our method for the simplest case, namely 
a fc* term in the dispersion relation. However, it looks 
like quite easy to include higher order derivatives. In 
order to do that, one should compute terms like ip^^P\ 



for p integer, by recursion and by exploiting the 0{n) 
invariance of the Green's functions, in the same way as 
we calculated Eq. In the case p = 3, we obtain 

few more terms which can be absorbed by a redefinition 
of the operator D. Therefore, the formal expression ((54|) 
does not change. A generic dispersion relation of the form 
!F{k'^), such that T is analytic, can always be expanded as 
a power series in . Therefore, it is reasonable to expect 
that Eq. ([54)l does not change, provided D is adequately 
defined. We hope to prove this conjecture in a future 
work. 
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APPENDIX A: EVALUATION OF I,{y,s) 
To evaluate leiVi s), we first write it as 

^^(y-^) = /(^e"'^'"^''^^""'"' 7 = -^^. (Al) 
By expanding the second exponential, we can write 

m 

The integral is now a product of n Gaussian integrals, 
hence 

(A3) 

By recalling the definition of Hermite polynomials 

^w = (-i^^'|r^"^'' (A4) 

and by changing variable, we finally find Eq. (j65p . 
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